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DURATION 
METHOD OF USE 



CAPACITY 



PUNCHING OF TAPES 



UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER 



LIBRARY ROUTINE K 5 - 182 



Autocorrelations 
Entire Program 
See below 

The program is read into the memory in the usual way * 
After the master tape is read in -a sum cheek is performed. 
If the master tape has "been read in incorrectly , ten 
sexadecimal characters will be punched. The master tape 
should then be read into the memory again, Next^ the 
data tape is placed in the reader and read in. The first 
character that appears on the print -out gives the type . 
of autocorrelation* This is followed by the mean of the 
observations and then the variance. Two columns appear 
next. The first column is the lag and the second is 
the autocorrelation. -After all the autocorrelations 
have been computed (the lags having increments of unity 
from first to last),, a new .pro^iem can be begun. 
The maximum .number of observations- this- program will 
handle is 795* 

Each observation which must be in the range *» 1 < x < 1 
is punched as a sign followed by up to 12 decimal digits. 
The final fraction in the sequence must be followed 
by an I, J F/ or L. Each character signifies a different 
type of autocorrelation' These are given by the 
following; 



N: 



1 v (x. * x) (x. , - x) 

~* L v i l+h 
n ..■ 



i=l 



with x. = x. if j > n, 
where s = 1 y (x. - xj 

5 i=i 



-2- 

1, ^ (x. - I) (x i+h - S 2 ) 

n-h ■' . ., 

s 

■j a "- n 

where x n = ~*"t* E x. ^ 
1 n-h . .. i ■ 
i=l 

and X,~ = ■ - -, Zj X . , » 

2 n-h . , l+h * 

i=l 

n n-h (x. - x) (x. - x) 

„ 1 v i l+n 

Fs — - L — — —77 

n-]a ... o 2_ 

i=l s 



x n-h (x i - ^Kx -= x 2 ) 

L: -^-r 2j ■<-* .-v v - ■-. ' ■ V. ' . :?! r-*p 



n~h 



1=1 s 1 .6 



— 2 

, n-h (x. -*, X.. ) 
where s, =* y v l . 1 

±=1 n-h 

— 2 
, n-h (x. ,, - x ) 
and s = y i+h .2 / 

i=l n-h 



This is followed "by an integer giving the initial lag 
desired^ a fifth hole character, an integer giving 
the last lag desired^ another fifth hole character ^ and 
finally by the number €>f decimal places to "be printed. 
The number of places to be printed must not exceed 9» 
When the problem is compietedj a new data may be read 
into the computer^ a different type of autocorrelation 
using the same data may be computed $ or the same type 
may be computed over a different range of lags» To 
use the data in the computer again it is only necessary 
to read in new parameters- 
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Time per Autocorrelation of lag a 

Type Computation 

N 2.8 n + 85.3 milliseconds 

J 4.l(n-h) +86.9 " 

F ^.l(n-h) + 87.2 M 

L 58.7(n-h) + 220 M 



Output 



50 (p+*0 milliseconds 



where n = number of observations 

and p = number of places each correlation is to be printed. 
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: 0EBER 
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i 



Routine X 1 

00 166K 
Roa&ine N 2 

00 192K 
Routine P 1 

00 220K 
Routine R 1 

00 3K 

00 F 

00 229F 

00 F 

00 166F 

00 F 

00 192F 

00 F 

00 220F 

00 25K 

92 135F 

41 5E 
50 S3 
50 1L 
26 S4 
40 3F 
L5 8s4 
10 20F 

42 13TL 
LO l40L 
40 16F 
LO 131L 
40 17F 
19 38F 
66 17F 

S5 S3 

40 18F 

L5 7L 
42 10L 

41 8F 



Decimal Order Input 



Input a Sequence of Decimal Fractions 



Print One Number 



Square Root Routine 



Read in Fractions 



Test how many numbers have been 
read in 



n 



l/n 



LOCATION 


ORDER 
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10 


4l 9F 








L5 F 


by 9, 16} from 17 




11 


40 IF 
50 IF 






12 


7J 18F 
40 2F 






13 


L4 8F 


< 






40 8F 




Compute mean and sum of squares 


l4 


50 2F 
7 J IF 






15 .. 


l4 9F 






t 


40 9F 






16 


F5 10L 
42 10L 






17 


LO 1371 








32 10L 




i 


18 


50 8f 
79 8f 






19 


L4 9F 








4o 9F 




Compute variance 


20 


L5 3F 
L4 24l 






21 


42 26L 
00 36F 






22 


82 4f 

1 




Punch out type of autocorrelation 




L5 3F 






23 


F4 231 








42 90L 




Determine type of autocorrelation 


24 


4l 13F 
4l 12F 




desired 


25 


4l 14F 
4l 15F 






26 


15 132L 






27 

j 


40 F 
15 1331 
4-2 311 


by 21 





i 

| XOOAIIGN 


ORDER 
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| 28 

i 
i 
i 

| 


50 5F 
81 ^F 


from 3^ 




29 

i 
1 


7^ 133L 
00 k-F 


f rom 30 




30 


91 ^F 








36 29L 




Read in Initial and last lag 


31 


S5 F 








40 F 


by 27, 32 




32 


F5 31L 
\2 31L 






33 


IxL JLi 

1*0 L 






3^ 


36 28L 

^1 6f 






35 


81 4f 

00 20F 




Head in number of places to be printed 


36 


k6 6f 
10 20F 






37 


IA 6f 
10 39F 






38 


IA- 6F 








-Jk 133L 




Compute print routine parameters 


39 


K5 F 






ko 


h6 A8L 
00 20F 






kl 


k6 118L 
92 135F 






k2 


F5 IIP 
IA 131L 






^3 


^0 11F 
92 515F 






hh 


92 987F 
L5 8F 






^ 


52 F 


"by 39 




* 


50 k5h 

i ; 




Print out mean 



LOCATION 


s 1 

ORDER 
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k6 


2.6 S5 
92 131F 






47 


92 987F 
L5 9F 






kQ 


52 F 


by 40 






50 48L 




Print out variance 


k9 


26 S5 








L5 134L 




1/10 


50 


4o 6f 








50 18F 


from 121 




51 


Fl 12F 
32 55L 






52 


L5 10F 
LO 17F 






53 


kO OF 
09 IF 






5^ 


10 39F 
66 OF 






55 


L5 10F 
kO 7F 






56 


S5F 








ko 19F 




Compute l/n-h or l/n 


57 


kl 20F 
4l 21F 






58 


kl 22F 
4l 23F 






59 


41 24F 








92 131F 






60 


L5 137L 
42 136L 






61 


92 515F 
L5 135L 






62 


42 63L 
L4 10F 






63 


42 65L 








L5 F 


by 62 , 8l| 








from 87 


i T 



, 1. 

LOCATION 


ORDER 


KOTES 
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6k 


10 IF 




V 2 








kO OF 








65 


32 65L 
L5 F 


by 63, 83 








66 


10 IF 
kO IF 




x i + h/ 2 






61 


50 19F 
7 J OF 










68 


kO 2F 
50 2F 










69 


7 J IF 
IX 21F 










70 


40 21F 




^ x * x i+h where 5 


= for r circular 






Fl 12F 




4p(n - Sh) 


= 1 for r mon*clreylar 




71 


32 80L 
L5 2F 










72 


IA 22F 
40 22F 




Ex. 

1 










2(h-hJ 






75 


50 19F 
7 J IF 










Jk 


40 3F 
I&;23F 










75 


40 23F 
L3 15F 




Zx i+h . / 2(n~h) 






76 


32 80L 
50 2F 




Zx 2 / 2(n-li) 






77 


7J OF 
IA 24F 










78 


kQ 24F 
50 3F 










79 


7 J IF 
iA 20F 




Zx 2 . / 2(n~h) 






80 


40 20F 
F5 63L 










81 


42 63L 












F5 65L 

i 






t 


- 


s 











LOCATION 
82 

83 

85 
86 
87 



89 
90 



LO 136L 
32 8*kL 
l4 136L 
il-2 65L 
22 85L 
L5 135L 
^2 65L 
F5 7F 
kQ 7F 
LO 17F 
32 63L 
L5 9F 
10 2F 
lj-0 TF 
L5 8F 
10 IF 
h-0 k-F 
2.6 F 



91 


50 k¥ 




22 102L 


92 


50 22F 




22 108L 


93 


L5 22F 




22 103L 


9H- 


50 22F 




79 22F 


95 


Lk 2^+F 




kO 2k¥. 


96 


50 23F 




79 23F 


97 


IA 20F 




kO 20F 


98 


50 20F 




75 2*KF 


99 


40 IF 




S5 F 



from 82 
from 84 



by 23; 
from 130 



from 108, 102 
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Are computations complete? 



Store s /k Ijsi J 



Store x/2 in h 



Betermine t^pe of eorrelatioii 





t 




LOCATION 


ORDER 
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100 


ho OF 
50 100L 






101 


22 S6 
ko 7F 






102 


26 92L 








79 ^F 


from 91 




103 


26 109L 








LO 4-F 


from 93 




104 


40 OF 
L5 23F 






105 


LO 1*-F 
kO IF 






106 


50 OF 
7 J IF 






107 


lh 21F 
40 21^ 






108 


26 92L 








79 23F 


from 92 




109 


lA 21F 
kO OF 


from .103 




110 


50 OF 
75 6F 






111 


kO OF 
L7 F 






"112 


12 JF 








36 127L 




Test for sealing 


113 


L5 OF 
66 7F 






Hk 


S5 IF 
^0 *KF 






U5 


32 H5L 








L5 10F 




Friiit out lag 


116 


50 kk~F 
50 U6L 






U7 


26 S5 









L5 ^F 







K 5 



LOCATION ORDER 

118 52 F ^ "by 4l 
50 118L 

119 26 S5 
F5 1QF 

120 40 10F 
LO 11F 

121 52 50L 
24 122L 

122 92 155F 
81 4F 

125 LO 158L 
40 5F 

124 56 20L 
L5 159 

125 4o 8s4 
L5 5F 

126 L4 25S4 
22 10S4 

127 50 6F from 112 
7J 154L 

128 40 6F 
L5 H8L 

129 L4 h4l 
46 U8L 

150 22 90L 
00 F 

151 80 F 
00 F 

152 LL 4095F 
LL 4095F 

155 00 F 

00 10F 

154 00 F 
00 1000 0000 OOOOJ 

155 00 S5 
00 S5 
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Print out autocorrelation 



Test if new data or new parameters 
have "been read 



LOCATION 



136 

137 
138 

139 
]AO 



ORDER 



32.65L 
L5 F 
SI 9F 
L5 F 
00. F 
00 12F 
'kO S3 
L5 8s4 
80 F 
40.S3 
,00 500K 
Routine X 7 
24 25N 
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by 4 



Sum dneck 



*fc- 



■*r 



